home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cgeev.z / cgeev
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCGGGGEEEEEEEEVVVV((((3333FFFF))))                                                            CCCCGGGGEEEEEEEEVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CGEEV - compute for an N-by-N complex nonsymmetric matrix A, the
  10.      eigenvalues and, optionally, the left and/or right eigenvectors
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR, LDVR, WORK,
  14.                        LWORK, RWORK, INFO )
  15.  
  16.          CHARACTER     JOBVL, JOBVR
  17.  
  18.          INTEGER       INFO, LDA, LDVL, LDVR, LWORK, N
  19.  
  20.          REAL          RWORK( * )
  21.  
  22.          COMPLEX       A( LDA, * ), VL( LDVL, * ), VR( LDVR, * ), W( * ),
  23.                        WORK( * )
  24.  
  25. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  26.      CGEEV computes for an N-by-N complex nonsymmetric matrix A, the
  27.      eigenvalues and, optionally, the left and/or right eigenvectors.
  28.  
  29.      The right eigenvector v(j) of A satisfies
  30.                       A * v(j) = lambda(j) * v(j)
  31.      where lambda(j) is its eigenvalue.
  32.      The left eigenvector u(j) of A satisfies
  33.                    u(j)**H * A = lambda(j) * u(j)**H
  34.      where u(j)**H denotes the conjugate transpose of u(j).
  35.  
  36.      The computed eigenvectors are normalized to have Euclidean norm equal to
  37.      1 and largest component real.
  38.  
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      JOBVL   (input) CHARACTER*1
  42.              = 'N': left eigenvectors of A are not computed;
  43.              = 'V': left eigenvectors of are computed.
  44.  
  45.      JOBVR   (input) CHARACTER*1
  46.              = 'N': right eigenvectors of A are not computed;
  47.              = 'V': right eigenvectors of A are computed.
  48.  
  49.      N       (input) INTEGER
  50.              The order of the matrix A. N >= 0.
  51.  
  52.      A       (input/output) COMPLEX array, dimension (LDA,N)
  53.              On entry, the N-by-N matrix A.  On exit, A has been overwritten.
  54.  
  55.      LDA     (input) INTEGER
  56.              The leading dimension of the array A.  LDA >= max(1,N).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCGGGGEEEEEEEEVVVV((((3333FFFF))))                                                            CCCCGGGGEEEEEEEEVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      W       (output) COMPLEX array, dimension (N)
  75.              W contains the computed eigenvalues.
  76.  
  77.      VL      (output) COMPLEX array, dimension (LDVL,N)
  78.              If JOBVL = 'V', the left eigenvectors u(j) are stored one after
  79.              another in the columns of VL, in the same order as their
  80.              eigenvalues.  If JOBVL = 'N', VL is not referenced.  u(j) =
  81.              VL(:,j), the j-th column of VL.
  82.  
  83.      LDVL    (input) INTEGER
  84.              The leading dimension of the array VL.  LDVL >= 1; if JOBVL =
  85.              'V', LDVL >= N.
  86.  
  87.      VR      (output) COMPLEX array, dimension (LDVR,N)
  88.              If JOBVR = 'V', the right eigenvectors v(j) are stored one after
  89.              another in the columns of VR, in the same order as their
  90.              eigenvalues.  If JOBVR = 'N', VR is not referenced.  v(j) =
  91.              VR(:,j), the j-th column of VR.
  92.  
  93.      LDVR    (input) INTEGER
  94.              The leading dimension of the array VR.  LDVR >= 1; if JOBVR =
  95.              'V', LDVR >= N.
  96.  
  97.      WORK    (workspace/output) COMPLEX array, dimension (LWORK)
  98.              On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
  99.  
  100.      LWORK   (input) INTEGER
  101.              The dimension of the array WORK.  LWORK >= max(1,2*N).  For good
  102.              performance, LWORK must generally be larger.
  103.  
  104.      RWORK   (workspace) REAL array, dimension (2*N)
  105.  
  106.      INFO    (output) INTEGER
  107.              = 0:  successful exit
  108.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  109.              > 0:  if INFO = i, the QR algorithm failed to compute all the
  110.              eigenvalues, and no eigenvectors have been computed; elements and
  111.              i+1:N of W contain eigenvalues which have converged.
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.